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Abstract: ANXA2{AxmQxm.A2), a calcium-dependent phospholipid bind¬ 
ing protein, is involved in various Ca 2+ -related biological activities. In the 
present study, full-length cDNA of ANXA2 was isolated from the velvet 
antler tip tissue of sika deer (Cervus nippon hortulorum); the amino acid 
sequence and gene expression was analyzed by using bioinformatics and 
real-time reverse transcriptase polymerase chain reaction (RT-PCR) 
techniques. Nucleotide sequence analysis reveals that the full-length 
cDNA of the ANXA2 gene was 1372 bp, of which 1020 bp was in the 
open-reading frame (ORF) encoding 339 amino acids; its relative mo¬ 
lecular weight was 38.3 kDa; and isoelectric point was 6.72. Sequence 
analysis indicates that the protein includes four conserved tan¬ 
dem-duplication ANX domains. The gene-accession nucleotide sequence 
number in GenBank is JX315571. Expression analysis by RT-PCR re¬ 
veals that ANXA2 gene expression has a significant positive correlation 
with the antler-tissue mineralization process, indicating that this gene 
may play an important role in the regulation of antler-tissue mineraliza¬ 
tion. 

Key words: velvet antler, AnnexinA2 (ANXA2) gene, cDNA library, 
clone, real-time quantitative RT-PCR 

Introduction 

Velvet antler is the dense fuzzed, immature horn of male deer 
(excluding reindeer) (Brockes and Kumar 2002; Price and Allen 
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2004), it has attracted the attentions of many scientists due to its 
incomparable growth rate, regeneration mechanism, and per¬ 
ceived magical medicinal value (Li et al. 1995; Odelberg 2005). 

Velvet antler has a long growth period of about 90-120 days 
every year, and the fastest growth rate has been recorded at 1-2 
cm per day (Brockes and Kumar 2005; Han et al. 2005). Some 
researchers divide the development of antler into two phases: 
growth period and ossification period (Gerke et al. 2005; 
Markoff and Gerke 2005). In the first phase, growth is dominant, 
with a slow ossification rate. In contrast, during in the second 
phase, antler can be ossified fast, being accompanied with the 
accumulation of large amounts of mineral matters; however, this 
process produces an inhibitory effect on antler growth. Antler 
with a high level of ossification contains a high concentration of 
mineral matter and a declining ratio of organic components. 
Therefore, ossification degree indexes the old and young degree 
of velvet antler and the changes of components, further decides 
the rank of velvet antler. Improvement in the yield of velvet ant¬ 
ler and the economic benefits of the deer industry should not 
only inhibit or reduce antler ossification but also not influence 
antler growth rate. Illuminating the relationship between the 
growth rate and ossification is the basis to achieve this objective 
(Benaud et al. 2004; Rescher and Gerke 2004). 

ANXA2 is a calcium/phospholipids-binding multi-functional 
protein (Filipenko et al. 2004; Yamada et al. 2005; Mussunoor 
and Murray 2008). It plays important roles in the formation of 
biomembranes, endocytosis/exocytosis, osteoblast and osteal 
reabsorption, membrane transport, ion channels, DNA synthesis, 
and cell proliferation (Chiang et al. 1999; Gerke and Moss 2002; 
Okuse et al. 2002). In the present study, the ANXA2 gene in the 
constructed velvet antler tip tissue of sika deer belongs to this 
high-abundance expression gene, indicating that it may be an 
important factor regulating antler development. 

To date, there is no record of deer-derived ANXA2 gene se¬ 
quences in GenBank. In this paper, the full-length cDNA of 
ANXA2 genes were cloned successfully from the full-length 
cDNA library of velvet antler tip tissues of sika deer. Further, the 
gene structure and the expression levels of different tissue layers 
of the antler tip were studied using bioinformatics and real-time 
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quantitative RT-PCR technique. These results provide an impor¬ 
tant basis for further study of the biological functions of ANXA2 
genes, molecular mechanism of regulating antler development, 
and velvet antler production. 

Materials and methods 

Velvet antler tissue 

Antler grown for 80 days were collected from anaesthetized 
adult sika deer (Cervus nippon hortulorum) stags. The distal 4 
cm of the tips was removed. Different tissue layers (skin, reserve 
mesenchyme, precartilage, and cartilage) of the growing antlers 
were determined as described by Li et al. (2002). After being 
harvested, all samples were immediately preserved in liquid 
nitrogen and kept at -80°C until they were used for isolating the 
RNA. 

Methods for obtaining full length cDNA of ANXA2 genes 

Total RNA extraction and cDNA library construction were based 
on the protocols of SV Total RNA Isolation System reagent 
(Promega) and Creator™SMART™cDNA Library Construction 
Kit (Clontech), respectively. Using Ml3 primer, single colonies 
randomly selected from velvet antler tip tissue in the original 
library were put through large-scale 5’-end EST sequencing, by 
using a 3730XL sequencer (Applied Biosystems, Foster City, 
CA, USA). High-quality ESTs were clustered and spliced, using 
Phrap software. The consensus sequence splice was programmed 
by BLASTX and BLASTN; from there, the sequence homology 
comparison and function notes of the obtained genes were con¬ 
ducted with non-redundant proteins and nucleic acid database in 
GenBank. ANXA2 genes with corresponding positive colonies 
were selected to carry out PCR identification, bacterial amplifi¬ 
cation, and plasmid extraction, and then two-way sequencing for 
the plasmid DNA. 

Nucleotide sequence analysis 

The ANXA2 gene sequence was analyzed and compared using 
the BLAST P and ORF search programs with GenBank database 
search. The multiple sequence alignment of ANXA2 gene was 
carried out, using the Clustal W analysis program; the sig¬ 
nal-peptide site was predicted by Signal P3.0; and the ANXA2 
protein MW and pi were computed by using the ProtParam tool. 
A phylogenetic tree, based on evolutionary distances, was con¬ 
structed from amino-acid sequences using the njplotWIN95 pro¬ 
gram. 

Quantification of the ANXA2 gene expression by real-time 
RT-PCR 

Total RNAs using SV Total RNA Isolation System reagent 
(Promega) were isolated from different tissues including skin, 
mesenchyme, precartilage, and cartilage. The DNA residue was 
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removed by DNase I digesting, at 37°C for 30 min. Four micro¬ 
gram of the total RNA were used in each lane and electrophore- 
sed in a 0.8% agarose gel, at 100 V/12 cm for 15 min. 
First-strand cDNA synthesis was perfonned using M-MLV re¬ 
verse transcriptase (TaKaRa Biotechnology (Dalian) Co., Ltd. 
Japan) to transcribe poly (A) + RNA with oligo-d(T)18 and ran¬ 
dom six as the primers; reaction conditions were set as the 
manufacturer’s instructions. The cDNA was used for the assay of 
quantitative real-time PCR. The SYBR Green I real-time 
RT-PCR assay was carried out in an Option-II Sequence Detec¬ 
tion System (MJ Research, U.S.). 

In a real-time RT-PCR study, specific primers (ANXA2- F: 
5' - AGCTGC AGG A A AT CAACAG-3' and ANXA2- R: 

5'-ATCAATGACAGAGCCATCC-3') were used to amplify a 
146 bp fragment with cDNA from skin, reserve mesenchyme, 
precartilage and cartilage, and organs using 18S as a positive 
control. Quantitative real-time PCR primers were designed on 
the basis of EST sequences of the ANXA2 gene and 18S rRNA 
gene from deer-antler tip tissue by using Primer 5.0 and Oligo 
6.0 software. The EST sequences of the two genes were obtained 
by large-scale EST sequencing of the full-length cDNA library 
from deer-antler tip tissue constructed by our lab. The amplifica¬ 
tions were perfonned in a 96-well plate, in a 25 pL reaction 
volume containing 12.5 pL of 2 x SYBR Green Master Mix 
(TARAKA), 2.5 pL (each) of ANXA2-F and ANX42-R primers 
(10 mM), 1 pL of template, and 9 pL of DEPC-water. The ther¬ 
mal profile for SYBR Green real-time PCR was 95°C for 2 min, 
followed by 45 cycles of 94°C for 12 s, and 55°C for 30 s. In the 
96-well plate, each sample was amplified in triplicate. 
DEPC-water for the replacement of template was used as nega¬ 
tive control. The relative expression was calculated as 2 AACt 
(Livak and Schmittgen 2001). 

Statistical analysis 

A multiple comparisons (Duncan’s) test was conducted to com¬ 
pare significant differences in the ANXA2 expression between 
skin, mesenchyme, precartilage, and cartilage, using the 
SPSS 13.0 software. A significance level of p = 0.05 was chosen. 

Results and discussion 

Sequencing and bioinformatics analysis of ANXA2 

Computer analysis, using the BLAST algorithm, confirmed that 
the selected sequence corresponded to ANXA2. The full-length 
ANXA2 cDNA contains 1372 bp, of which 21 bp in the 
5'-untranslated region (UTR); 1020 bp in the open reading frame 
(ORF); and 331 bp in 3'-UTR with poly (A) tail (Fig. 1). The 
ORF encodes a polypeptide of 339 amino acids. The calculated 
molecular mass of the mature protein (339 amino acids) is 38.3 
kDa, with an estimated pi of 6.72. ANXA2 is in shortage of signal 
peptides, so they belong to non-secretory proteins and can not be 
excreted out cells by the classical protein excretion pathway. 
However, in a temperature stress state, these substances can be 
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transferred from the cytoplasm to the membrane surface (Deora NCBI GenBank as accession number JX315571. 
et al. 2004). The cDNA sequence has been submitted to the 


1 aattcggccattacggccgggsfgtctaccgttcatgaaattctgtgcaagctcagtttg 60 
MSTVHETLCKLSL 
61 gagggtgatcactccacacctccaagtgcatacgggtcagtcaaagcatacactaacttt 120 
EGDHSTPPSAYGSVKAYTNF 
121 gatgctgagcgggatgctctgaacattgaaacagccatcaagaccaaaggtgtggatgag 180 
DAERDALNTETAIKTKGVDE 
181 gtcaccatcgtcaacatcctgaccaaccgcagcaatgagcagagacaggatattgccttc 240 
VTIVNILTNRSNEQRQDIAF 
241 gcctaccagaggaggaccaagaaggaacttgcatcagcactgaagtcagccttgtccggc 300 
AYQRRTKKELASALKSALSG 
301 cacctggagacagtgattttgggcctattgaaaacacctgctcagtatgatgcttccgag 360 
HLETVILGLLKTPAQYDASE 
361 ctgaaagcgtccatgaaggggctagggactgatgaggactctctcattgagatcatctgc 420 
LKASMKGLGTDEDSLIEI IC 
421 tcaaggaccaaccaggagctgcaggaaatcaacagagtctacaaggaaatgtacaagacc 480 
SRTNQELQEINRVYKEMYKT 
481 gatctggagaaggacatcgtttccgacacatctggcgacttccgcaagctgatggtcgcc 540 
DLEKDIVSDTSGDFRKLMVA 
541 ctcgcgaagggtcggagagcagaggatggctctgtcattgattatgaactgattgaccag 600 
LAKGRRAEDGSVIDYELTDQ 
601 gacgccagggatctctatgatgctggtgtgaagaggaaaggaactgatgttcccaagtgg 660 
DARDLYDAGVKRKGTDVPKW 
661 atcagcatcatgaccgagcggagcgtgtgccacctccagaaagtatttgaaaggtacaag 720 
ISTMTERSVCHLQKVFERYK 
721 agctacagcccttatgacatgctggagagcatcaagaaggaggtcaaaggagacctggaa 780 
SYSPYDMLESIKKEVKGDLE 
781 aatgcctttctgaatctggtccagtgcattcagaacaagcccctgtattttgctgacaga 840 
NAFLNLVQCIQNKPLYFADR 
841 ctgtatgactccatgaagggcaagggcactcgcgacaaggtgttgattagaatcatggtc 900 
LYDSMKGKGTRDKVLIRIMV 
901 tcgcgcagtgaagtggacatgttgaaaattagatctgaattcaagaagaagtacgggaag 960 
SRSEVDMLKTRSEFKKKYGK 
961 tccctctactactacattcagcaagacaccaagggcgactaccagaaagcgctgctgtac 1020 
SLYYYIQQDTKGDYQKALLY 
1021 ctgtgtggtggggatgacfgaacccctgcgacggcccgagcgtccaggacggctgctgcg 1080 
L C G G D D * 

1081 tgctccagctaccggttctccacaatcagcttgcagctaacagccccctgtgcccatccc 1140 
1141 tgggaggatgacattagcgctgcccaccccgcccgctcttttgttagtttagtttcctac 1200 
1201 gcatfacetggctttcccggctagtccctctttatagccaaagaaatgaacattccaagg 1260 
1261 attt ggaagtg aaatctatgatgtgaaacactttgcctctcgagtactgtgtcataaacg 1320 
1321 ggtgaataaactgaatttttaccaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 1372 


Fig. 1 : Nucleotide and deduced amino acid sequences of ANXA2 cDNA of velvet antler from sika deer 


The C'-end of ANXA2 gene has four homologous duplication bound with Ca 2+ , phospholipids, and F-actin, thus having impor- 
sequences (ANX domains) that are from the amino acids 37 to tant significance in cytoskeleton recombination and various 

102, 122 to 174, 207 to 259, and 282 to 334 (Fig. 2). They can be membrane-related changes (Hayes et al. 2006). 


Query seq. 

Non-specific 

hits 

Superfanilies 



Fig. 2: Search for the conserved domains in deduced amino acid sequence of ANXA2 cDNA 


The N'-end of the gene contains three loci, namely Tyr23, PDGF receptors. Therefore, the three protein kinases play im- 

Serll, and Ser25. Among them, Seri and Ser25 can be phos- portant roles in cellular signal transduction and cell differentia- 

phorylated by the action of protein kinase; Tyr23 can be modi- tion (Babbin et al. 2007). 

fied by some membrane-binding kinases such as the IGF-1 and 
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Homology comparison of ANXA2 

The comparison of the ORFs with other known ANXA2s indi¬ 
cated that the ANXA2 shows a homology of Identities = 338/339 
(99%) with Bos Taurus, Identities = 333/339 (98%) with Homo 
sapiens, Identities = 323/339 (95%) with Rattus sp., Identities = 
304/339 (89%) with Callus gallus. Identities =275/339 (81%) 
with Xenopus laevis. Identities =215/339 (63%) with Monop- 
terus albus. Homological analysis indicated that ANXA2 gene 
encoding regions of different genus show a high homology. 
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Multiple alignment of ANXA2 protein sequences showed large 
variation in the beginning of C'-end, but in general, ANXA2 ma¬ 
ture protein is highly conservative (Fig. 3). The evolutionary tree 
can be divided into two branches: one branch is mammals; the 
other branch contains birds, amphibians, and fish. The mammal 
graph showed that in the gene locus, sika deer has the closest 
evolutionary relationship with cattle (Bos Taurus) and the fur¬ 
thest genetic relationship with fish such as Monopterus albus 
(Fig. 4). This result is in agreement with their traditional taxo¬ 
nomic position for the tested species, according to animal evolu¬ 
tion relationships. 


★ • - *•* •*•* _ ^ ^ * -*• • k k ************ 

MSTVHEILCKLSLEGD—HSTPPSAYGSVKAYTNFDAERDALNIETAIKTKGVDEVTIVN 58 

MSTVHEILCKLSLEGD-HST PPSAYGSVKAYTNFDAERDALNIETAIKTKGVDEVTIVN 58 

MSTVHEILCKLSLEGD-HSTPPSAYGSVKAYTNFDAERDALNIETAIKTKGVDEVTIVN 58 

MSTVHEILCKLSLEGDSQHSTPPSAYGSVKPYTNFDAERDALNIETAIKTKGVDEVTIVN 60 

MSTVHEILSKLSLEGD-HSLPPSAYATVKAYSNFDADRDAAALEAAIKTKGVDEVTIIN 58 

MALIHEILGKLSLEGN—QSSSRQSKLGSVKAATHFDAEKDAAAIETAIKTKGVDELTIIN 59 
MALVSEFLSQLTLAHT GEEEPKFPTWPVRDFDPAKDAARIEAAIKTKGVDEQTIID 57 

k k k ^ k k k k • k k k k • * • ★ ★ • ★ *k»k.kkkk k k • • • k k k • * • kkkkkkkkk»k 

ILTNRSNEQRQDIAFAYQRRTKKELASALKSALSGHLETVILGLLKTPAQYDASELKASM 118 
ILTNRSNEQRQDIAFAYQRRTKKELASALKSALSGHLETVILGLLKTPAQYDASELKASM 118 
ILTNRSNAQRQDIAFAYQRRTKKELASALKSALSGHLETVILGLLKTPAQYDASELKASM 118 
ILTNRSNAQRQDIAFAYQRRTKKELPSAMKSALSGHLETVMLGLFKTPAQYDASELKASM 120 
ILTNRSNEQRQDIAFAYQRRTKKELSAALKSALSGHLEAVILGLLKTPSQYDASELKAAM 118 
ILTNRSNEQRQDIAFAFHRRTKKDLPSALKGALSGNLETVMLGLIKTRPQYDASELKASM 119 
ILTRRSYEQRRDIAFEYERLAKKDLITALKGALSGSLEALILGLMKSTAQYDASELKASM 117 

★ k kk»kkk*kk»kkkkkkkkkk''k .*;;;**; kkkkkk»kkkkkkkk 

RAEDGSviDYELIDQDARDLYDAGVKRKGTDVPKWISIMTERSVCHLQKVFERYKSYSPY 238 
RAEDGSVIDYELIDQDARDLYDAGVKRKGTDVPKWISIMTERSVCHLQKVFERYKSYSPY 238 
RAEDGSVIDYELIDQDARDLYDAGVKRKGTDVPKWISIMTERSVPHLQKVFDRYKSYSPY 238 
RAEDGSVIDYELIDQDARELYDAGVKRKGTDVPKWISIMTERSVCHLQKVFERYKSYSPY 240 
RCEDTSVIDYELIDQDARELYDAGVKRKGTDVPKWINIMTERSVPHLQKVFERYKSYSPY 238 
RQEDGNMVDYEKIDQDARELYEAGVKRKGTDVTKWITIMTERSISHLQKVFERYKSYSPY 239 

RDEPSNWDYQKIDDDARS LYEAGVKRKGTDVATWNSIFSQRSIPHLQKVFERYKSYSPY 237 

★ - ★**•********. •**_**•*» ■ * * • • ** ^ ** • • - ★ - . ^ ★ - 

DMLESIKKEVKGDLENAFLNLVQCIQNKPLYFADRLYDSMKGKGTRDKVLIRIMVSRSEV 298 
DMLESIKKEVKGDLENAFLNLVQCIQNKPLYFADRLYDSMKGKGTRDKVLIRIMVSRSEV 298 
DMLESIRKEVKGDLENAFLNLVQCIQNKPLYFADRLYDSMKGKGTRDKVLIRIMVSRSEV 298 
DMLESIRKEVKGDLENAFLNLVQCIQNKPLYFADRLYDSMKGKGTRDKVLIRIMVSRSEV 300 
DMLESIKKEVKGDLENAFLNLVQCIQNKQLYFADRLYDSMKGKGTRDKVLIRIMVSRCEV 298 
DIEESIKKEVKGDLENAFLNLVQCIQNKPLYFADRLYESMKGKGTKDKILIRIMVSRCEL 299 
DMKESIRKEVKGDLEKSFLTLVECFENKQLYFANRLNEAMKSKGAKEKWTRIMVSRCEV 297 

k • • k k • kkk • • • ★ ★ • k • kkkkkk»kkk k k k k • k 

DMLKIRSEFKKKYGKSLYYYIQQDTKGDYQKALLYLCGGDD 339 
DMLKIRSEFKKKYGKSLYYYIQQDTKGDYQKALLYLCGGDD 339 
DMLKIRSEFKRKYGKSLYYYIQQDTKGDYQKALLYLCGGDD 339 
DMLKIRSEFKRKYGKSLYYFIQQYTKGDYQKALLYLCGGDD 341 
DMLKIKSEFKRKYGKSLYYFIQQDTKGDYQRALLNLCGGED 339 
DMLKIRQEFKKKYGKSLHYFIGQDTKGDYQRALLNLCGGDD 340 
DLMKIRTEFKRHQKRSLYQAIAEHTKGDYQKALLSLCGGDD 338 


Fig. 3: Multiple alignment of ANXA2 protein sequences 


Tissue expression of ANXA2 

ANXA2 expressed in each tissue of velvet antler tip is shown in 
Fig. 5. The real-time RT-PCR showed that the ANXA2 gene was 
detected in skin, reserve mesenchyme, precartilage, and cartilage. 
However, the gene has a low expression level in the skin layer 
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and reserve mesenchyme layer of velvet antler, but a high ex¬ 
pression level in the precartilage and cartilage layers. The gene 
expression has a significant positive correlation with antler tissue 
mineralization process, indicating an important role of the gene 
in the process of antler tissue mineralization. 

ANXA2 genes belong to protein super-family members that 
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can be combined to negative-charge membrane phospholipids by 
Ca 2+ mediated pathway, thus these proteins are involved in vari¬ 
ous Ca 2+ -related biological activities. Francis and Suttie (1998) 
detected ANXA2 gene expression in red deer antler by using the 
Western blotting technique, and found that ANXA2 proteins were 
located in the edge of osteoblastic membrane. Based on our 
findings and previous results, it is assumed that ANXA2 genes 
play a partial role in antler-tissue mineralization by participating 
in the formation of a Ca 2+ -ion channel. 


Monopterus albus 

Fig. 4: The phylogenetic tree of ANXA2 from animals. Clustal W was 
used to establish the phylogenetic tree, and the result was displayed using 
Treeview software. 

5 
4 
3 
2 
1 
0 

12 3 4 

Different tissues of velvet antler tip 

Fig. 5: Expression of ANXA2 in different tissues of velvet antler tip, 
1-skin; 2-reserve mesenchyme; 3-precartilage; 4-cartilage 

Velvet antler, as the most active growing point, stores large 
numbers of cell factors; these substances are the biochemical 
basis for clinical functions of velvet antler. Currently, our labo¬ 
ratory is conducting studies on the biological functions of 
ANXA2 genes. Research on antler ossification can not only reveal 
the old and young degree of velvet antler and component 
changes, which directly decide the optimal picking time and rank, 
but it also has significance for the construction of antler bone 
development models. Moreover, the related researches are help¬ 
ful to improve velvet antler yield and quality, and provide the 
theoretical foundation for new drug development. 
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